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Fusion rules in turbulence address the asymptotic prop- 
erties of many-point correlation functions when some of the 
coordinates are very close to each other. Here we put to exper- 
imental test some non-trivial consequences of the fusion rules 
for scalar correlations in turbulence. To this aim we examine 
passive turbulent advection as well as convective turbulence. 
Adding one assumption to the fusion rules one obtains a pre- 
diction for universal conditional statistics of gradient fields. 
We examine the conditional average of the scalar dissipation 
field (V 2 T(r)|T(r + R) - T(r)) for R in the inertial range, 
and find that it is linear in T(r + R) — T(r) with a fully de- 
termined proportionality constant. The implications of these 
findings for the general scaling theory of scalar turbulence are 
discussed. 
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The aim of this Letter is to present new analysis of ex- 
perimental data JIJ-pl pertaining to turbulent scalar ad- 
vection and to discuss the implications of this analysis in 
the context of fusion rules and conditional averages. We 
begin with a short theoretical background to the issues 
in order to make this Letter self contained. Turbulent 
advection is described mathematically by the equation 
of motion for a scalar field T(r,t) 



[fit + u(r , t) ■ V] T(r, t) = «V 2 T(r , t) , 



(1) 



where k is the scalar diffusivity and u(r,t) is the turbu- 
lent velocity field responsible for the advection of T(r, t). 
The problem of "passive" scalar advection is the one in 
which the properties of u(r, t) are not affected by those 
of the scalar T(r,t). In "active" scalar problems, like 
turbulent convection, the velocity field and its statisti- 
cal properties are coupled with those of the scalar field, 
and Eq. (|l|) has to be supplemented with an additional 
equation for u(r,t) and T(r,t). In our thinking below 
we consider passive as well as active scalar fields. In both 
cases we are interested in the limit of large Peclet num- 
ber Pe which is defined as UlL/k where Ul is the typical 
velocity difference across the outer scale L of turbulence. 

The statistical properties of the scalar fields are com- 
monly discussed in terms of the so called "structure func- 
tions" S 2n (R) defined as 



S 2n (R) = ([T(r + R,t)-T(r,t)r n ) 



(2) 



where (. . .) stands for an ensemble average. In writing 
this equation, we assume that the statistics of the veloc- 
ity field leads to a stationary and spatially homogeneous 



ensemble of the scalar T. If the statistics is also isotropic, 
then S 2n becomes a function of R only, independent of 
the direction of R. The scaling exponents of the struc- 
ture functions S 2n (R) characterize their R dependence in 
the limit of large Pe, 



S 2n (R)ocR^, 



(3) 



when R is in the "inertial" interval of scales that will 
be discussed later in this paper. One of the fundamental 
questions in the theory of turbulent advection is what the 
numerical values of the exponents Q 2n are, and whether 
they conform with classical Kolmogorov type arguments, 
or rather exhibit the phenomenon of multiscaling. 

An important equation to analyze in this context is 
the so-called "balance equation" which is obtained by 
writing Eq. (Q) twice at points r and r + R, subtracting 
the equations, and multiplying the result by 2n[T(r + 
R,t) — T(r, i)] 2 " -1 . Taking the ensemble average and 
using the symmetry between the two points analyzed, 
one finds the balance equation 



D 2n {R) = J 2 n{R) , 



(4) 



where D 2n (R) stems from the convective term in ([j]) and 
Jin{R) stems from the diffusion term: 

J 2n (R) = -Ann (V 2 T(r)[T(r + R) - T(r)] 2 "- 1 ) . (5) 

It was argued recently, first by Kraichnan jj] and later 
in refs. that balance equations play a very impor- 

tant role in providing non-perturbative relations that can 
determine, or severely constrain, the values of the scaling 
exponents C, 2n - A good example is Kraichnan's model of 
passive scalar advection || , in which the velocity field u 
is delta-correlated in time, but exhibits power law scal- 
ing in space. In this case the convective term Din can be 
calculated exactly in terms of Sin j§, 

D 2n (R) = -R^- d -^R d ^h{R)^Ls 2 n{R) , (6) 

with d being the space dimension and h(R) the scalar 
part of the eddy-diffusivity, h(R) cx R^ h with £/, a scal- 
ing exponent. If we could represent exactly also the RHS 
J 2n (R) in terms of S 2n (R), we could evaluate all the scal- 
ing exponents Q 2n from the balance equation (Q). Here 
is where the fusion rules come in. The fusion rules ap- 
pear naturally in the analytic theory of Navier-Stokes 
turbulence|7]|^-|llj and passive-scalar turbulent advec- 
tion Bpd|,|12|, and they determine the analytic structure 
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of n-point correlation functions when a group of coor- 
dinates tend towards each other. In the case of scalar 
advection we consider simultaneous many-point correla- 
tion functions of field differences: 

^2n{r a \ri,r 2 . ■ .r 2n ) = (7) 
(6T(r , n) ST(r , r 2 ) . . . 6T{r , r„)> . 

We note that the previously defined structure functions 
S 2n are obtained by "fusing" all the coordinates n ■ ■ ■ T% n 
to one coordinate Tq + R. The fusion rules were derived 
in p"T[ for systems that enjoy universality of the scaling 
exponents (i.e the scaling exponents do not depend on 
the detailed form of the driving of the turbulent flows) , 
and whose correlation functions T 2n are homogeneous 
functions of their arguments, 

^2n(Ar |Ari, . . . , Xr 2n ) = \ C2n J 7 2n (r \r 1 , . . . ,r 2 „) . (8) 

This form applies whenever all the distances \r$ — ro| 
are in the so-called "inertial range", between the outer 
scale L and the appropriate dissipative scale of the sys- 
tem, denoted below as r\. The fusion rules address the 
asymptotic properties of T 2n when a group of p points, 
p < 2n — 1 tend towards r (\ri — r \ ~ p for all i < p), 
while all the other coordinates remain at a larger dis- 
tance R from ro (\ri — Vq\ ~ R for i > p, and R 3> p). 
In particular under the two general assumptions of scale 
invariance and universality of the scaling exponents the 
fusion rules state that to leading order in p/R 

^ 2 „(r \r + p,r + R,..., r + R) ~ ^j^S 2n (R). 

D2{-tt) 

(9) 

This forms holds as long as p is in the inertial range. 

We show now how to use this fusion rule to calculate 
J2n(R)- First write it as 

J2n{R) = -4k« lim V 2 JF 2 n(ro\r +p, r a +R, r +R). 

In using the fusion rule (^) to evaluate this quantity we 
interpret the limit p —> as a limit p — > rj. This seems 
natural for large Peclet numbers when 77 — > 0. It is im- 
portant however to stress that there is a hidden assump- 
tion here. We expect the function J : 2 n { r o\fa + P- r o + 
R, . . . , ro + R), which is a function of p and R to change 
its analytic behavior as a function of p. This change oc- 
curs at the viscous crossover scale rj. The issue is whether 
this crossover scale is n and R independent. That this 
is so has been proven for Kraichnan's model of turbulent 
advection (ll]] but not in general. We believe that this is 
more generally true due to the linearity of the equation of 
motion ([|), independently of the statistical properties of 
the driving velocity field. The experimental results which 
we discuss later in this Letter will strongly indicate that 



this is the case in a wide context of scalar turbulent fields. 
We caution the reader that this is not so in Navier-Stokes 
turbulence. With this in mind we write 

hn{R) ~ -^n[V 2 p S 2 {p)\ p ^}S 2n {R)/S 2 {R) . (10) 

Using the fact that the mean of the scalar dissipation 
field, denoted e, is evaluated as I ~ n[V p \S 2 (p)\ p = v ] , and 
also the fact that in the inertial range J 2 (R) = — 4e, we 
write 

J 2n (R) = nC 2n J 2 S 2n (R)/S 2 (R) , (11) 

where C 2n is an as yet unknown dimensionless coefficient, 
butC 2 = 1. Eq.(|ll|) was suggested for Kraichnan's model 
in p| and derived in 0. Here we propose that it holds 
in a much wider context. To this end we turn now to the 
analyzes of experimental data. 

We first display experimental results that confirm the 
theoretical prediction (|TT|). The results show that to a 
good accuracy C 2n w 1 for all n and R. The theoretical 
consequences of this n and i?-independence of C 2n will 
be discussed after examining the data. 

FIG. 1. A plot of log|J an (fl)/(2nK)| vs. 

log I (2k) _1 J2S 2 n(R)/S2(R) \ for n = 2 (squares), 3 (triangles), 
4 (diamonds), 5 (stars), and 6 (circles) and R in the inertial 
range. The data are taken from Yale [1]. The line is not a fit, 
but the theoretical expectation with slope 1 and intercept 0. 

We use temperature data measured in the wake of a 
heated cylinder jij. Water of speed 5 m/s flowed past 
a heated cylinder of diameter 19 mm (Reynolds number 
= 9.5 x 10 4 ). The temperature was measured at a fixed 
point downstream of the cylinder on the wake center- 
line. The cylinder was heated so slightly that the buoy- 
ancy term was unimportant and temperature acted as 
a passive scalar. Temperature was measured as a func- 
tion of time, and we use here the standard Taylor hy- 
pothesis that surrogates time derivatives for space deriva- 
tives. In Fig. 1 we display J 2n (R) / (2nn) as a function 
of (2k) -1 J 2 S 2n (R) / S 2 (R) for n varying from 2 to 6, and 
for various R values in the inertial range. We see that 
all the points fall on a line whose slope is unity to high 
accuracy, and whose intercept (in log-log plot) is very 
closely zero. This good agreement is a confirmation of 
the validity of the fusion rules. In addition, this agree- 
ment lends support to the assumption that rj is n and R 
independent. It should be stressed that individual tests 
at various values of n as a function of R corroborate the 
same conclusion, i.e. Eq.(p"i"|) is supported by the ex- 
perimental data with C 2n being near unity. The most 
sensitive test of the alleged constancy of the coefficients 
C 2n is obtained by dividing J 2n (R) by nJ 2 S 2n (R) / S 2 (R) 
for all the available values n and R. The result of such a 
test is shown in Fig. 2. We see that all the measured val- 
ues of C 2n are concentrated within the interval (0.75, 1) 
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for all separation within the inertial interval. Consider- 
ing the fact that the quantities themselves vary in this 
region over five orders of magnitude, we interpret this as 
a good indication for the independence of C 2n of R and 
n. The R independence is very clear, and is a direct test 
of the fusion rules. The weak n dependence seems to in- 
dicate that Can decreases slightly with n; this may arise 
from the limited accuracy of the data. We are reluctant 
to make a strong claim about the accuracy of 10'th or 
12'th order structure functions. 

FIG. 2. A detailed test of the coefficient Cin, see text for 
details. The symbols are the same as in Fig.l. The small 
systematic decrease of Cin with n may be due to insufficient 
accuracy at the tails of the probability distribution which be- 
come more important at large values of n. 

Let us accept for now the evidence that the coefficients 
Gi n in Eq. (|TT|) arc ri-independent, and look for a way to 
understand it. Note that J 2n {R) can be exactly written 
in terms of conditional averages in the form 

J 2n (R) = -4nn J dST(r, r + R)P[ST{r, r + R)] (12) 

x {5T(r,r + R)} 2 "- 1 (\7 2 T(r)\6T(r,r + R)) , 

where 5T(r,r + R) = T(r + R) - T(r). We see that 
the conditional average (\7 2 T(r)\ST(r, r + i?)]} appears 
as a natural object that needs to be determined. If we 
make the assumption that the conditional average, which 
in general is a function of the two variables R and ST is 
factorizable as a function of R times a function of ST, 
then the only possible such form is 

-4«(V 2 T(r)|<5T(r,r + #)) = -J*6T(r,r+R). (13) 

With this form in @ we regain the RHS of Q, but 
since the conditional average cannot be a function of n 
the coefficient C 2n mustbe n- independent. We note that 
Kraichnan conjectured that the conditional average is 
linear in ST(r,r + R) in the context of the Kraichnan 
model, and numerical simulations supporting this con- 
jecture were presented J^]. Moreover, linearity approx- 
imations on conditional average of the form (S7 2 X\X) 
and relation like ([ll]) have also been studied earlier by 
Ching jl3| and by Pope and Ching P4| . In this Letter 
we propose that the linearity of the conditional average 
in ST(r,r + R) is a general property of a wider variety 
of turbulent advection problems. From our discussion, it 
is clear that other results on conditional statistics can be 
derived in a similar fashion. 

FIG. 3. The conditional average in Eq.(16) as measured 
from the Yale data [1] normalized by the measured value of 
Jz/Sz^K) as a function of 8T(r, r + R) for three different val- 
ues of R. The different R values are designated by triangles, 
squares and circles respectively. 



Before we proceed to the implications of ( |13[ ) we 
present the experimental evidence for its validity. In 
Fig. 3 we present results from the same data set that was 
used above. We show the conditional average as a func- 
tion of 5T(r,r + R) for various values of R. The line 
passing through the data points is not a fit, but rather 
the line required by Eq.(|l3|). We note that points be- 
longing to different values of R fall on the same line, 
indicating that indeed the conditional average is a func- 
tion of ST(r, r + R) times a function of R, and that we 
identified correctly the function of R as J2/S 2 (R). To 
test the generality of this result we analyzed a second 
data set from the convective hard turbulence regime of 
the well documented Chicago experiment j^P]. The ex- 
periment was performed in a cylindrical box of helium 
gas heated from below, and the Rayleigh number can be 
as high as 10 15 . The box has a diameter of 20 cm and 
a height of 40 cm. The temperature at the center of the 
box was measured as a function of time, and we use the 
same Taylor-hypothesis to analyze the conditional aver- 
age. The results are shown in Fig. 4. Although we see 
larger statistical scatter at the ends of the plot, the basic 
assertion of linearity with the correct slope are confirmed. 

FIG. 4. Same as Fig. 3 but computed from the Chicago 
data, Refs. 2,3. 

As explained, the linearity of the conditional average 
in ST, Eq.(p^|) was not derived from first principles. To 
stress the theoretical interest in such a derivation we con- 
sider briefly another form of J2n{R) that is obtained by 
moving around one of the gradients in (||) . Up to a term 
that is negligible for R in the inertial range we can write 

J2n(R) = -4n(2n - l)n (\X7T(r)\ 2 T(r + R) - T(r)} 2 "- 2 ) 

(14) 

Accepting Eq.([y]) with C 2n = 1 we can write 

- AK{\VT{r)\ 2 [T(r + R) -T{r)] 2n - 2 ) 

= J 2 S 2n (R)/(2n-l)S 2 (R). (15) 

The LHS can be written, similarly to (|l2|), in terms of 
the conditional average (| VT(r)| 2 |5T(r, r + R)). It is 
obvious however that now we cannot assume that this 
quantity factorizes into a function of ST times a function 
of R. If it did, the dependence on ST must have been 
(ST) 2 in order to give us S 2n (R) on the RHS of @. But 
we can never obtain in this way the explicit l/(2n — 1) 
factor. This underlines the fact that the factorization in 
( |l3| ) is far from being obvious or trivial. Currently we do 
not know the deep reason why the conditional averages of 
V 2 T afford factorization. We pose this as an important 
issue for further theoretical research. 

Lastly, we comment on the implications of these find- 
ings for the exponents . As discussed above, if we know 
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the functional form of Jin and the coefficient, we can use [15] K. Gawedzki and A. Kupiainen, Phys. Rev. Lett. 75, 

the balance equation (||) to compute the scaling expo- 3608 (1995). 

nents, provided that we know the nonlinear term Dm ■ In 

the context of the Kraichnan model the latter is known 

exactly, and the balance equation leads to a quadratic 

equation for the exponents C„, with the solution 

(2n = kl(2-d+ ^{C2+dY+±dUn-l)\ ■ (16) 

These are the exponents that were conjectured by Kraich- 
nan. We note that these exponents are in disagreement 
with the calculations of refs. fl^,[l^] which attempted to 
compute the exponents by perturbative methods using 
as a small parameter either 2 — £2 or the inverse dimen- 
sion 1/d. If it turns out indeed that (|l6| ) is the correct 
non-perturbative result, one needs to carefully rethink 
the meaning of these perturbative calculations. 
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